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Abstract 

We study the effect of intense magnetic fields on the phase diagram of cold, strongly interacting 
matter within an extended version of the Nambu-Jona-Lasinio model that includes flavor mixing 
effects and vector interactions. Different values of the relevant model parameters in acceptable 
ranges are considered. Charge neutrality and beta equilibrium effects, which are specially relevant 
to the study of compact stars, are also taken into account. In this case the behavior of leptons is 
discussed. 
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I. INTRODUCTION 


The influence of intense magnetic fields on the properties of strongly interacting matter 
has become an issue of increasing interest in recent years Jl]. This is mostly motivated by 
the realization that in some relevant physical situations, like high energy non-central heavy 
ion collisions j^j] and compact stellar objects called magnetars [3], very strong magnetic 
fields may be produced. Since in these systems extreme temperatures and/or densities may 
be found, it is interesting to investigate which modifications are induced by the presence of 
strong magnetic fields on the whole QCD phase diagram. Unfortunately, even in the absence 
of those fields, the present knowledge of such phase diagram is only schematic due to the 
well-known difficulty given bv the so-called sign problem which affects lattice calculations 

4], Of course, the presence of strong magnetic fields makes the 
situation even more complex. Thus, most of our present knowledge of their effect comes from 

and refs, therein). 


at finite chemical potential 


investigations performed in the framework of effective models (see e.g. jfj] 

In this contribution we present some results of a study of the phase diagram of cold quark 
matter subject to intense magnetic fields in the framework of a generalized Nambu-Jona- 
Lasinio (NJL) model. The NJL-type models are effective relativistic quark models for non 


out and interactions are 
6] it only includes scalar 


perturbative QCD, where gluon degrees of freedom are integratec 
modelled through point like interactions. In its simplest version 
and pseudo scalar interactions that describe chiral symmetry breaking effectively. As well 
known, however, a more detailed description of the low-energy quark dynamics requires that 
other channels like flavor mixing and vector meson interactions are taken into account [/J. 
In fact, some aspects of the effect of those interactions on the magnetized quark matter have 


already been investigated 


8, isl. The purpose of the present work is to extend those analyzes 
by performing a detailed study of the resulting cold matter phase diagrams, including their 
dependence on the parameters that regulate the strength of these interactions. Moreover, 
the behavior of cold magnetized quark matter under conditions relevant for the physics of 
compact stars will also be considered. One of the phenomena to be discussed in detail 


is that related to the so-called inverse magnetic catalysis 


temperature and moderate values of the chemical potentials 10]. It is important to mention 


MC) expected to exist at low 


that, despite bearing a similar name, this phenomenon is different from the inverse magnetic 
catalysis at finite temperature found in lattice QCD (LQCD). Concerning the latter one, we 
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recall that most effective models foresee that at zero chemical potential a crossover transition 
is obtained at a pseudo critical temperature T c that increases with an increasing magnetic 
field, a behavior which is contrary to the one found in LQCD calculations 11]. In fact, 
recently there have been significant efforts to modify the models such that they in corp orate 
a mechanism that could lead to inverse magnetic catalysis around T c (see Ref. 12] for a 
recent review on this issue). We should stress, however, that this is not expected to affect 
the low temperature behavior discussed in this work. 

This paper is organized as follows. In Sec. II we provide some details of the model and 
its parametrizations as well as the way to deal with an external constant magnetic field. In 
Sec. Ill we present and discuss our results for symmetric quark matter. The situation for 
stellar matter is analyzed in Sec. IV. Finally, our conclusions are given in Sec. V. 


II. FORMALISM 

We consider a generalized NJL-type SU(2) Lagrangian density which includes a scalar- 
pseudoscalar interaction, vector-axial vector and the t'Hooft determinant interaction [7]. In 
the presence of an external magnetic field and chemical potential it reads: 


£. = 0 (ip — m c + 0 7 0 ) 0 + £, 


int 


(i) 


where 


int 


6 6 

Y [{^ T ^Y + O/’WsTaV’) 2 + G 2 Y [('07/VaV’) 2 + (07^0) ‘ 


a =0 


a =0 


+G 3 [(^ 7 ^) 2 + (07/G50) 2 + G 4 [(07^0) 2 - (07m7s0) 2 + 2 G D (, d + + cL) (2) 


T 

Here, G* with i — 1,4 and Go are coupling constants, -0 = {u,d) represents a quark field 
with two flavors, d± = det [0(1 ± 75 ) 0 ], 0 = diag {fi U} fid) the quark chemical potentials, 
m c is the (current) mass matrix that we take to be the same for both flavors, 77 = /, where 
/ is the unit matrix in the two flavor space, and r a , 0 < a < 3 denote the Pauli matrices. 
The coupling of the quarks to the electromagnetic field is implemented through the 
covariant derivative D M — iqA^ where q represents the quark electric charge matrix 

q = diag (q u ,qd) where q u / 2 = — qj = e/3. In the present work we consider a static and 
constant magnetic field in the 3-direction, A^ = d^XiB. In the mean-field approximation 
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the associated grand-canonical thermodynamical potential for cold and dense quark matter 
reads 


qMfa = _Y^6 f + Q pot (3) 

}=u,d 

where Of gives the contribution from the gas of quasi-particles of each flavor / = u,d and 
can be written as the sum of 3 contributions [lsj] 


N 

QVSLC _ 1 ' C 


f 8 tt 2 
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A (2A 2 + M )) ^ A 2 + M 2 - Mj In 
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amed c 

/ 4tj-2 


C ( 1 , 0 ) (-l,a;/) - ~{x) - x f ) In x f + ^ 


..max 

v f 


u=0 


MB E Ijlfjfl 2 - s f(u, B ) 


-Sf(u, B) 2 In 


0f + y/ 2 / - s /0> 5 ) 5 

s f (v,B) 


(4) 


where Mj = m c + crf and fij = fif — uif, with oq- and (U/ being the mean held values of the 
scalar and vector meson holds, respectively. A represents a non covariant ultraviolet cutoff 
and ('( 1 ’ 0 )(—1 , Xf) = d((z,Xf)/dz\ z= -i where ((z,Xf) is the Riemann-Hurwitz zeta function. 
In addition, Sf(u,B ) = \JM 2 + 2\qf\Bv while Xf = Mj/(2\qf\B). In QJ ed , the sum is over 
the Landau levels (LLs), represented by u, while a u = 2 — 5 u0 is a degeneracy factor and 
u m ax ^j ie igj-ggg^ integer that satishes < (Jij — M 2 )/(2\qf\B). 

The Q pot contribution reads 

(1 - C s )(cr 2 + CT 2 ) - 2c s a u a d (1 - 2c„)(u; 2 + uj 2 d ) + c v uj u uj d 


Qpot 


(5) 


% s (l-2c s ) 8g v (1 — 2c v ) 

where we have introduced a convenient parametrization of the coupling constants in terms 
of the quantities g s , c s , g v and c v given by 


9s 

c s 


G i + Gd 


G D 

Gi + G D 


9v — G2 + G3 + G4 

_ G3 + G 4 
2 {Cj 2 + G3 + G4) 


The relevant gap equations are given by 


dQ MFA (cr u ,cr d ,UJ u ,UJ d ) 
d(a u ,a d ,uj u ,uj d ) 


( 6 ) 

(7) 
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The solutions to the gap equations minimize the thermodynamic potential with respect 
to the quark masses Mf , but the ojf derivatives amount to a consistency condition and 
the potential is actually a maximum with respect to these variables. Several solutions will 
generally exist, corresponding to different possible phases, and the most stable solution is 
that which minimizes the thermodynamic potential with respect to Mf. 

In our calculations we will consider first the simpler case of symmetric matter where 
both quarks carry the same chemical potential //. Afterwards, we will analyze the case of 
stellar matter in which leptons are also present and /3-equilibrium and charge neutrality are 
imposed. In this case the chemical potential for each quark, /if. is a function of quark number 
chemical potential /i and the lepton chemical potentials which have to be self-consistently 
determined. 

In order to analyze the dependence of the results on the model parameters, we will 
consider two SU(2) NJL model parameterizations. Set 1 corresponds to that leading to 
M 0 = 340 MeV while Set 2 to that leading to M 0 = 400 MeV. Here, M 0 represents the 
vacuum quark effective mass in the absence of external magnetic fields. The corresponding 
model parameters are listed in Table HI 


TABLE I. Parameter sets for the NJL SU(2) model. 


Parameter set 

M 0 

m 

fhA 2 

A 

- <uu > x / 3 


MeV 

MeV 


MeV 

MeV 

Set 1 

340 

5.595 

2.212 

620.9 

244.3 

Set 2 

400 

5.833 

2.440 

587.9 

240.9 


The presence of the t’Hooft determinant interaction is very important since it reflects 
the Ha( 1)—anomaly of QCD. Its strength, and consequently the amount of flavor mixing 
induced by this term, is controlled by the parameter c s . An estimate for its value can be 
obtained from the r] — rf mass splitting within the 3—flavor NJL model [m 1]. This leads 


to c, ~ 0.2 


15| . In any case, to obtain a full understanding of the effects of flavor mixing 


we will vary the value of c s in a range going from 0, which corresponds to a situation in 
which the two flavors are completely decoupled, to 0.5, being this the case of maximum 


flavor mixing described for example in Ref. 


16]. Regarding the vector coupling term it is 
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important to recall that one can obtain naturally the terms proportional to G 2 if one starts 
from a QCD-inspired color current-current interaction and then performs a Fierz transform 
into color-singlet channels, and that in this case the relation between coupling strengths is 
G 2 = G i/2 [7]. Yet, the value of g v cannot be accurately determined from experiments nor 
from lattice QCD simulations and this is why it has been taken as a free parameter in most 
works. In the present work we take 0 < g v /g s < 0.5. It is worth mentioning that due to the 
mixing of pseudoscalar and longitudinal axial vector interaction terms, pseudoscalar meson 
properties depend on G 2 . Thus, strictly speaking the parameters given in Table 1 only lead 


to the empirical values of f n and when G 2 = 0. However, as shown in Ref. 


13], ra* 


and f 7 r only change by ~ 10% when G 0 /G 1 increases from 0 to 1. Thus, for simplicity, we 
will keep the model parameter values fixed when varying g v . A last comment regarding c„, 
i.e. the parameter that regulates the ratio between the singlet and octet vector-axial vector 
interaction strengths: we will take it as a free parameter in the range 0 < c v < 1/2. Note 
that for c v = 1/2 only singlet vector-axial interactions are present and, thus, there is no 
mixing between the pseudoscalar and longitudinal axial vector channels. 

We end this section by describing the way in which the different phases of the magnetized 
quark matter will be denoted as well as the procedure used to identify the boundaries between 
them. For the phases we adopt the notation of Refs. flbl . nj. Thus, the vacuum (i.e. fully 
chirally broken) phase is denoted by B, the massive phases in which Mf depends on the 
chemical potential by C Q and, finally, the chirally restored phases by A a . Here, cc is a set 
of two numbers indicating the highest LL populated for each flavor. To obtain the critical 
chemical potentials at a given eB we proceed as follows. In the case of first order phase 
transitions we calculate the thermodynamical potential for each of the neighboring phases 
(that is, the global minimum with respect to Mf) and then search for the chemical potential 
values at which they become degenerate. In the case of crossover transitions the critical 
value is identified by the peak of the chiral susceptibility corresponding to each quark flavor, 
defined as d< 'ip'ijj >/dmf. 


III. RESULTS FOR SYMMETRIC MATTER 

In this section we present the results obtained for the case of symmetric matter. These 
results were obtained solving the set of coupled “gap equations” (1T4|) for different values of 
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magnetic field and chemical potential. 


A. Effect of the flavor mixing interactions 

To neglect vector interactions implies taking G 2 = G 3 = G 4 = 0 in Eq. (J2J), while G 1 ^ 0 
and Gd 7 ^ 0. Therefore, Eqs. (JTJ) , will become a set of two coupled equations that must 
be solved for the independent variables M u and Ma■ The parameter c s acts as a coupling 
between both flavors and we study how the phase transitions are modified as we vary this 
parameter in the range 0 < c s < 0.5. The case c s = 0.5 corresponds to ordinary NJL where 
flavor mixing is maximum and, thus, both flavors have identical behavior. In fact, the first 
term in Eq. (J5]) will tend to infinity as c s goes to 0.5 unless M u = Md , which leaves only 
one equation to be solved. But if c s < 0.5, then both masses will be independent variables, 
and transitions for each flavor might occur simultaneously or not in different regions of the 
phase diagrams. 

The phase diagrams for both parameter sets and several values of flavor mixing can be 
seen in Fig. [0 We will start by commenting some general features which are common to all 
phase diagrams discussed in this work. It is seen that chiral symmetry is completely broken 
for chemical potentials well below Mq and that restoration occurs for high enough chemical 
potentials, usually accompanied by a large drop in the dressed mass. The inclusion of a 
constant magnetic held modifies the quark dispersion relation, introducing Landau levels 
(LL’s) into its spectrum. A consequence of this is that chiral symmetry restoration might 
occur in several steps as chemical potential is increased, each of which is a transition where 
quark population appears on previously unoccupied LL’s. The restored chiral symmetry 
region consists of several phases with different number of LL’s occupied, which are separated 
by the so-called Van Alphen De Haas transitions, whose form is in the absence of vector 
mesons approximately /i c = ^ 2 k\qf\B (being this condition exact for the chiral case). As 
a result of the different quark electric charges, one up transitions is found every two down 
transitions when a phase diagram is traversed in the magnetic held direction at high enough 
fixed / 1 . In all phase diagrams, a “main transition” is found, which separates the vacuum 
phase from the phases with populated LL’s and at fixed eB it is the one with the lowest 
possible chemical potential (e.g.: lower black line in bottom left diagram in Fig. [[]). In some 
cases, another main transition within the populated phase exists. It connects phases with 
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partially restored symmetry and low LL population to the fully restored symmetry phases, 
where the Van-Alphen De Haas transitions are present and LL level population may be much 
higher. All of this gives rise to a potentially complex phase diagram, whose precise form 
depends on the parameter set and magnetic held (The B = 0 case is simple yet depending on 
the parameter set there can be a few differences). The lower main transition usually shows, 


for moderate magnetic he 


magnetic anticatalysis 


ds, a decrease of /i c when eB increases that is sometimes called 
101 ] (even though the name is more generally used to refer to the 
decrease of the quark condensates as eB increases). For higher values of the magnetic held 
this tendency is reverted, and this gives rise to a characteristic curve in the main transition 
line to which we will refer as the “IMC well”. It is also worth noting that some of the other 
transitions in the phase diagrams also exhibit an IMC-like behavior. 

Now we will discuss how the phase diagrams are modified as the parameter c s is varied. 
The upper two panels correspond to the c s = 0 case (no havor mixing), which was previously 
studied in Ref. jljij]. In a sense, each havor will have its own independent phase diagram 
because the gap equations are decoupled. Here, each line corresponds to a transition where 
LL population of a single quark havor occurs (red lines for down quarks, blue lines for 
up quarks). Both havors will coincide for eB = 0 where SU(2) symmetry is recovered 
and behave differently as eB increases, due to their different electric charges. Since this is 
the only difference between the equations for both havors and since it only appears in the 
product qjB, the down havor phase diagram may be obtained from the other one through the 
replacement q u B = (2 qd)B = q^B) which amounts to a stretching of the up havor phase 
diagram along the eB axis. In other words, since the down quark has a smaller coupling 
to the held than the up quark, it will require a magnetic held twice as large to replicate 
the effect on an up quark. As a consequence of this, the IMC wells are shifted with respect 
to each other, so there will be a large well-distinguished region where up quarks exist in 
the lowest LL (LLL) while down quarks are in vacuum (OB), and another region where the 
opposite occurs (BO). 

For Unite c s (second row onwards in Fig. |TJ) , the coupling between havors creates a complex 
pattern, where transitions move closer together to the point of coalescing in some regions, 
that is, the LL population changes simultaneously for both havors (these are represented by 
black lines). We can see that already for c s = 0.03, the transitions in both parameter sets 
occur together for low magnetic helds, and then separate for eB = 0.02 GeV 2 in Set 1 and 






eB = 0.1 GeV 2 in Set 2. Note that for c s = 0, the transitions for both flavors cross around 
eB = 0.12 GeV 2 (for both sets). When flavor mixing is introduced, this crossing point 
transforms into a line, that is, both transitions occur together once again in an interval of 
magnetic held values. When c s = 0.03, they separate at eB = 0.4 GeV 2 , which is intuitive 
since a higher magnetic held will further break SU(2) havor symmetry. For c s = 0.2 this 
separation is no longer seen in the diagrams but it can be guessed that it does occur beyond 
eB = 1 GeV 2 , and that for any value of c s < 0.5 there will always be a large enough magnetic 
held that will cause the u and d main transitions to separate. 

For a better understanding of the physical meaning of the transition lines, we present 
in Fig. [2] the dressed masses for both havors for Set 2, c s = 0.03, for eB = 0.11 GeV 2 . 
When the hrst discontinuity is encountered, at ji — 370 MeV, M u jumps to half its value 
and its LLL becomes populated. On the other hand, the down havor remains in vacuum. 
Actually, its mass presents a small discontinuity caused by the weak coupling to the up 
quark which is not to be interpreted as a down transition. The down quark LLL is occupied 
at n — 374 MeV. This is precisely the kind of behavior that generates a rich phase diagram 
for low c s values. The difference in masses is understood in terms of the effect of magnetic 
catalysis. Since the havors have different charges, they couple with different intensities to 
the magnetic held, so the up quark will have a larger mass in the vacuum phase and a lower 
one in the populated phases, which is consistent with the fact that mass increases with 
magnetic held in the vacuum phase and decreases in the populated phase. 

The behavior of the crossovers as c s is varied from 0 to 0.5 is interesting to note. For 
c s = 0, in Set 1, there is one crossover for each havor. When c s is increased, the up crossover 
is slightly shifted to the left when crossing from the 00 to the 01 phase, acquiring a small 
discontinuity. In turn, a down crossover starts to appear from the left in the latter phase. 
As c s is increased, these two crossovers move towards each other, while the up crossover in 
00 shifts to the right, also moving towards the down crossover originally existing in that 
same phase. The c s = 0.5 limit, in which down and up crossovers have joined, is achieved 
very slowly, being the crossovers still separated for c s ~ 0.4995. 

The hrst order transitions for both havors already occur together in the whole studied 
region for c s = 0.2, so the qualitative behavior is very similar to the full havor mixing case 
c s = 0.5. In fact, the model tends to full mixing quite quickly and only for c s < 0.1 are 
relevant mixture effects (or more precisely, the absence of it) actually seen. VA-dH lines are 
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brought together for a particularly small amount of mixing, already coinciding for c s = 0.03, 
while crossover transitions, on the other hand, tend much more slowly to the c s = 0.5 
behavior. 


B. Effect of the vector interactions 


In this section we analyze the effect of vector interaction terms. As discussed in the 
previous section, for c s > 0.1 the phase diagrams do not present qualitative variations. 
Here, therefore, we consider c s = 0.2 which also is in the range of realistic values suggested 
in Ref. 15]. Although the results to be shown below correspond to c v = 0, our studies show 
that only small quantitative differences occur when varying this parameter from 0 to 1/2. 

In Fig. [3]we present a series of phase diagrams obtained for different values of the ratio 
g v /g s - For Set 1 we observe that as g v /g s increases the two main transitions separate and 
several new transitions appear in between, in the low eB region of the diagram. These are 
partially restored symmetry regions, where quark mass acquires an intermediate value. For 
Set 2 there is a unique main transition for g v = 0, but already for g v /g s — 0.1 it will split 
into two, leaving a 00 phase in between which was not present before. For larger values of 
g v /g s the behavior is similar to that of Set 1. The existence of new transitions at low eB as 
g v /g s increases can also be appreciated in the left panel of Fig. [4j where we show d quark 
density normalized to nuclear matter density (po = 0.17 fm~ 3 ) as a function of chemical 
potential for eB = 0.016 GeV in Set 1. In the absence of vector interaction, the density 
jumps from close to 0 to 3 times nuclear matter density, while quark population jumps from 
00 to a phase where several LL’s are occupied. As g v increases, the amount of transitions 
increases too. In fact, increasing the vector interaction coupling has the same effect as going 
to a parameter set that reproduces a lower value of current mass Mo, described in Ref. 16]. 


Notice that the phase diagram for Set 2 and g v /g s = 0.3 is in this sense very similar to Set 1 
and g v /g s = 0. 

It is interesting to analyze the effect of the vector interactions on the so-called inverse 
magnetic catalysis (IMG) mentioned in the Introduction. We recall that the IMG is usually 
related to a decrease of the critical chemical potential at intermediate values of the magnetic 
fields, a phenomenon that can be clearly observed in all the phase diagrams plotted in Figs. [1] 
and [3] However, while from Fig. |T]we see that the variation of the strength of flavor mixing 


10 




interactions has basically no effect on the IMC effect, the situation is different for the vector 
interactions. In fact, from Fig. [3] we note that if we measure the depth of the IMC well as 
the difference in /x between the lowest critical chemical potential at vanishing magnetic field 
and the lowest possible one in the whole diagram, we find that this difference is reduced 
by an 84% for Set 1 and by a 67% for Set 2 when going from g v — 0 to g v /g s — 0.5. To 
explain this feature we recall that the IMC effect can be understood in terms of the extra 
cost in free energy to form a fermion-antifermion condensate at finite /x, Vt ext [loj]. In the 
absence of vector interactions this extra cost basically originates in the LLL contribution to 
the medium term. Considering the chiral limit for simplicity, this contribution can be shown 
to be proportional to B/j 2 for symmetric matter, and it tends to decrease the difference in 
free energy between the vacuum phase and the finite density phase. As it is clear from 
Eqs. dU [5]) the presence of the vector interactions introduces some modifications in Sl ext . In 
the case of the medium term, they imply the replacement fif —> jdf = fif — ujf. Moreover, 
there is a new contribution coming from £l pot . Thus, assuming as above that there is LLL 
dominance and that quarks are massless in the chirally restored phase, the extra cost for 
symmetric matter is 

^ ext = 4^2 + 4 — ( 8 ) 

where, for simplicity, we have assumed c v = 1/2. The generalization for arbitrary values 
of c v is straightforward and, in addition, the numerical dependence on c v of the estimates 
to be given below turns out to be negligible. Of course, uj should satisfy the associated 
gap equation which follows from Eq. (TT4T) in the Appendix. Within the above mentioned 
approximations, the solution of this equation is 


u = 


1 + 


7T 


N c gvY.f kt\ B 


-i 




( 9 ) 


Replacing in Eq. (JSJ) we get that Q ext can be expressed in a form similar to that ob¬ 
tained in the absence of vector interactions, but where /x has to multiplied by a factor 
l/y/l + N c g v eB/ tt 2 . Note that the actual values of the u — and d —quark charge have been 
already used to obtain this factor. Thus, in the region eB ~ 0.2 GeV 2 around which this 
expression is approximately valid we expect 


gc(.9v) 

MO) 



( 10 ) 
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For example, for Set 2 we obtain that the lowest critical chemical potential (which occurs 
at about eB = 0.24 GeV 2 ) increases by about 10 % when going from g v = 0 to g v /g s = 0.5. 
Together with the fact that the lowest p c at vanishing magnetic held stays basically constant 
for values of g v /g s >0.1 this explains the strong reduction of the IMC well. 

Another aspect of the reduction of the IMC phenomenon induced by the presence of 
the vector interactions can be observed in the right panel of Fig. |4j There, we plot the 
current quark mass for d quarks as a function of the magnetic held for Set 1, p, — 345 MeV 
and several g v /g s values, where the system is in the 00 phase at eB = 0. As discussed in 
Ref. Q i n this case an actual decrease of the mass as eB increases is expected to exist. As 
we see, however, such a decrease is slower for larger values of g v . It should be noticed that 
discontinuities appearing for g v /g s = 0 corresponds to the d quark transition from the phase 
00 to 01 and back. 


IV. RESULTS FOR STELLAR MATTER 


We now turn our attention to stellar matter, that is, matter where /3-equilibrium and 
charge neutrality are imposed. In this case, electrons and muons are introduced into the 
system so that the thermodynamical potential receives an extra contribution 13] 


S'" = EE 

l=e,n u=0 


\qi\B a v 
An 2 


Hi 


vV 2 - s iiy, B) 2 - si{y , B) 2 In 


Hi + \Zh\ ~ si(v,B¥ 
siiy,B) 


( 11 ) 


where = Int[(/i 2 — rri 2 )/(2\qi\B)) and Si(u, B) = \Jm 2 + 2\qi\Bu. We take m e = 

0.511 MeV and = 105.66 MeV. 

The ^-equilibrium and charge neutrality conditions read 


Hd Hu He j He Hfe 


( 12 ) 


and 

Pe + P[1 — 2 ( 2 p u — Pd) ) ( 13 ) 

respectively. The lepton densities appearing in the last equation can be easily obtained from 
the derivatives of the total thermodynamical potential with respect to the corresponding 
chemical potentials. 
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Following the discussions in the previous section only results for c s = 0.2 will be presented. 
If charge neutrality is imposed on our system, there will be a fixed relation between the 
densities of the up and down quarks, which will be necessarily different unless they are 
both zero. So, even though the value of c s = 0.2 was close enough to the full mixing case 
according to what was established in previous sections, charge neutrality will cause the 
flavors to behave differently among themselves. 

The phase diagrams in the /i — B plane for both parameter sets and for increasing values 
of g v are plotted in Fig. [5j The chemical potential g on the horizontal axis is now the 
quark number chemical potential, in terms of which the flavor chemical potentials read 
/i/ = g — g//i e when Eqs. (TT2|) are used. The introduction of stellar matter conditions has a 
few effects similar to those of vector interaction, in the sense that diagrams become similar 
to the ones corresponding to lower M 0 sets: main transitions separate for low eB and 
several transitions appear in the region between them. The magnetic anticatalysis effect 
is reduced also: The depth of the anticatalysis well, as defined in the previous section, 
is reduced from 35 MeV to 9 MeV in Set 1 and from 43 MeV to 22 MeV in Set 2. As 
discussed in Ref. Q this can be understood by generalizing to stellar matter the discussion 
given in Sec.IIIb. We see that again the minimum critical chemical potential occurs at 
eB ~ 0.2 GeV 2 . Around that value, and in the absence of vector interactions, the extra 
cost in free energy to form a fermion-antifermion pair at finite g is in this case proportional 
to Bp 2 with p 2 = l^/l/i 2 + /r 2 /3. Here, a generally small muonic contribution has been 

neglected. Using the relations obtained from the f3 equilibrium conditions, fif = g — qffi e , we 
get p = g(l — 2x/3 + 2x 2 /3) 1 / 2 , where x = /i e //i. Note that the minus sign in the (dominant) 
linear term follows from the fact that \q u \ = 2\qd\. The relevant value of x follows from the 
neutrality condition Eq. (TTH]1 . Assuming as before that in the chirally restored phase we are 
dealing with massless quarks one obtains x ~ 0.38 for /i ~ 350 MeV. Using this result we 
get p ~ 0.92 g. This implies that the extra cost in free energy is smaller than that required 
in the symmetric matter case for the same value of eB and g. Consequently, for a given eB, 
one needs a larger value of the chemical potential to induce the phase transition. In fact, we 
have /hf / li s c ym — 1.09 a value which is in good agreement with our full numerical results. In 
principle to determine the change in the IMC well we should also estimate the modification 
of ii c at eB = 0. As shown in Ref. [uj this value also increases when stellar conditions are 
imposed. However, such an increase is several times smaller than the one at eB 0.2 GeV 2 
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leading to a quenching of the IMC effect. As discussed in Sec.lllb, the introduction of vector 
interaction will further enhance these effects. Consequently we observe that anticatalysis 
has completely disappeared for g v = 0.5. 

The VA-dH transitions for different flavors also acquire relatively independent behaviors 
when charge neutrality is introduced. In particular, some down and up transitions occur 
simultaneously as we follow them upwards along the phase diagram and then separate visibly 
at a given chemical potential. This is most clearly seen at intermediate eB (~ 0.1 GeV 2 ) for 
the transitions separating the 02, 12 and 13 phases. The relationship for critical chemical 
potentials /i c = sj2 k\qj\B , which occurred for symmetric matter now holds for each flavor 
chemical potential separately, where these two are different between themselves and related 
through Eq. (fT2lh As a result of this, we see that in Fig. 0 which is instead plotted in terms 
of quark chemical potential, we roughly have one up transition every four down transitions, 
whereas in symmetric matter there was one up transition every two down transitions. In 
particular, for high g v , the transitions have a tendency to clump up in groups of three (two 
down transitions and one up) with an intermediate down transition which is well separated 
from this group. 

The transitions corresponding to the population of lepton Landau Levels are also included 
in Fig. [5j In all cases, it was seen that electron’s transition from vacuum to LLL occurs in 
the lower main transition, that is, quarks and electrons occupy their LLL simultaneously. 
This is understandable in that electrons have a very small, almost negligible mass, hence 
they will populate their LLL as soon as the associated chemical potential becomes finite. 
On the other hand, muon behavior is more complex. Since they have a larger mass than 
electrons it is expected that they will require a larger chemical potential to acquire a finite 
population. This is what actually happens in all our diagrams for low magnetic fields, 
where values larger than 430 MeV are needed for this, hence not appearing in the plotted 
region. However, at intermediate values magnetic held values (0.1 < eB < 0.2 GeV 2 ), 
this transition drops suddenly until it joins the lower main transition of the quarks. This 
happens because /i e increases with magnetic held, allowing for muon population at lower 
quark chemical potential. For higher eB values, all leptons and quark transitions occur 
together. The occupation of the hrst Landau Level of either lepton species requires much 
higher chemical potentials. This happens because leptons have a larger charge than quarks, 
leading to Landau Levels with larger energy, and because /i e is only a slowly increasing 
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function of p. The value needed for the first electron LL to be occupied, which is around 
250 MeV, is therefore only achieved near p = 900 MeV. In some phase diagrams, the muon 
transition is seen to affect the behavior of the crossovers. This occurs because the /i e is 
coupled to the order parameters, so a muon transition, with its corresponding jump in p e , 
will affect the values of the quark masses and its associated susceptibilities. For g v /g s = 0.5 
in Set 1, for example, the down crossover separating the 00 and 00 phases is smeared out near 
the muon transition. In other diagrams, like the one for g v /g s = 0.3, the muon transition 
has completely absorbed the quark crossover. 

In Fig. [6] we plot the normalized lepton densities together with those of the quarks as 
functions of the quark number chemical potential, for eB = 0.02 GeV 2 , g v /g s = 0.3 and pa¬ 
rameters of Set 1. All three densities become non zero simultaneously at p = 340 MeV, where 
they all occupy their corresponding LLL. However, while quark density grows monotonously 
after the first transition, it is seen that when the next quark transition is encountered p e 
decreases, and jumps to an even lower value in the transition that follows. Recalling that 
these correspond to transitions of d quarks and that the charge neutrality condition imposed 
is p e + Pn = (2p u — pd)/ 3, it can be understood that the electron density decreases whenever 
there is a d-transition and subsequent growth of pd (recall that = 0 in this range of chem¬ 
ical potential), without relevant changes in p u . This is the case of this range of chemical 
potential and this value of magnetic field. 

In Fig. [7] we show quark and lepton densities as functions of the magnetic field for p = 
0.36 GeV, for different strengths of the vector coupling. Results correspond again to Set 1. 
At low magnetic field, discontinuities in the quark densities, associated with the VA-dH 
transitions, are always encountered. For the cases with g v = 0 and 0.3 a region is found in 
which pd changes its monotony within a C-type phase. This is because quarks can be found 
in this phase for magnetic fields high enough so the mass has started to increase with eB, 
and is associated to the region of the phase diagram above the IMC well. This behavior is 
not observed for g v /g s = 0.5 where the IMC effect has disappeared. For eB ~ 0.2 GeV 2 , d 
quarks are already in a phase with a low k and u quarks with k = 0, all densities increase 
and it is interesting to notice how slopes change whenever a d transition is encountered. 
When d quarks finally reach the 0 phase all densities increase until they drop to zero when 
the last transition to vacuum is encountered. 
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V. CONCLUSIONS 


We have investigated how the phase diagram of cold strongly interacting matter is mod¬ 
ified in the presence of intense magnetic fields in the context of NJL-type models which 
include flavor mixing and vector interaction. The whole range of possible flavor mixing 
values was swept through, going from the situation in which the two flavors are completely 
decoupled, to the one in which they are fully mixed. For low mixing values, a complex phase 
diagram is generated, but already for c s ~ 0.1, phase diagrams display a behavior that is 
qualitatively very similar to the full mixing case. Since SU(3) estimates suggest a value of 
c s ~ 0.2, it can be concluded that the realistic flavor mixing range is similar in behavior to 
the full mixing case. In what followed, vector interactions and stellar matter conditions were 
introduced. The most notable effect observed is that they attenuate the Inverse Magnetic 
Catalysis phenomenon, to the point that it completely disappears when both effects are 
jointly taken into account. It was also found that introducing vector interaction causes the 
two main transitions to separate and additional phases to appear between these two, effects 
that are similar to those occurring when changing to a parameter set that fits to a smaller 
dressed mass. In the stellar matter case, the behavior of the leptons was also studied. It was 
found that while electron LLL becomes populated simultaneously with quarks, the muon 
transition presents a more complicated dependence with the magnetic held. Namely, muon 
LLL requires a very high chemical potential to become populated at low eB , however, it 
joins the main quark transition together with electrons for high enough magnetic held. 
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APPENDIX 


The explicit form of the gap equations is 
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FIG. 1. (Color online) Phase diagrams in the eB — [i plane for different values of flavor-mixing 
parameter c s . To simplify the figure we have introduced a compact notation to indicate the phases. 
The pair of integers mn corresponds to the C mn phase and the pair fhn to the A mn phase. The 
case in which one quark is in a C-type phase and the other in the A-type phase is indicated by 
putting a bar on top of the integer associated with the A-type phase. 
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FIG. 2. (Color online) Dressed masses for both flavors for Set 2, c s = 0.03, for eB = 0.11 GeV 2 . 
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FIG. 3. (Color Online) Phase diagrams in the eB — n plane for different values of g v /g s ■ Different 
phases are denoted as in Fig. |Tj 
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FIG. 4. (Color online) Quark density over nuclear matter density as function of the chemical 
potential for eB = 0.016 GeV 2 (left) and current mass as function of magnetic field for g = 
0.345 GeV (right) for different values of g v /g s ■ Results correspond to d flavor and were obtained 
with parameter Set 1. 
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FIG. 5. (Color online) Phase diagrams in the eB — g plane for stellar matter and different values of 
g v /g s . Different phases are denoted as in Fig. [I] The pink dotted line represents muon transition 
from vacuum to LLL. 
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FIG. 6. (Color online) Quark and electronic densities over nuclear matter density as functions of 
the chemical potential for eB = 0.02 GeV 2 . Results were obtained with g v /g s = 0.3 and parameter 
Set 1. 


24 









12 


Pu . 


Pd . / 


4 Pe - / 


' 4 Pu .- / 

- 

/ / 
/ / , 

- g v ig s = 0 / // 


y'x// 

■ 

-v" ... / 

■ 

__ . _ . _ 

■ 


4 

2 

0 



0.01 0.1 

eB [GeV 2 ] 


FIG. 7. (Color online) Quark and lepton densities over nuclear matter density as function of the 
magnetic field for g = 0.36 GeV for different values of g v /g s ■ Results were obtained with parameter 
Set 1. 
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